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Abstract 

A  study  has  been  carried  out  to  assess  the 
effect  of  flow  nonuniformities  in  AEDC  Tunnel  B  on 
surface  pressures  of  slender  bodies.  The  approach 
taken  was  to  use  flow  profiles  measured  in  a  recent 
calibration  program  as  inflow  boundary  conditions 
for  a  CFD  solution  to  compare  with  the  idealized 
case  where  a  uniform  free  stream  is  assumed.  The 
results  of  this  study  indicate  that  flow  nonuniformi¬ 
ties  are  at  least  partially  responsible  for  discrepan¬ 
cies  observed  at  low  Reynolds  number  (Re/L  =  1 .0 
x  106/ft)  on  slender  axisymmetric  bodies  at  zero 
angle  of  attack.  There  is  no  corresponding  effect  at 
high  tunnel  Reynolds  number  (Re/L  =  3.5  x  106/ft). 
At  high  Reynolds  number,  for  angles  of  attack  that 
displace  the  nose  of  the  vehicle  more  that  6-7  in. 
from  the  tunnel  centerline,  there  is  a  small  increase 
in  surface  pressure  caused  by  a  radial  variation  in 
Mach  number. 

Nomenclature 

Cp  .  Pressure  coefficient,  (P  -  P„)/(0.5pQM) 
e  -  Internal  energy/unit  volume,  f^/sec2 

Re  -  Reynolds  number,  dimensionless 

L  -  Body  length,  in. 

M  -  Mach  number,  dimensionless 

P  -  Pressure,  psi 

Q  -  Total  velocity,  ft/sec 

S  -  Circumferential  distance,  in. 


T 

-  Temperature,  'R 

U 

-  X-  component  of  velocity,  ft/sec 

V 

-  Y-component  of  velocity,  ft/sec 

w 

-  Z-component  of  velocity,  ft/sec 

X 

-  Axial  Model  Coordinate,  in. 

Y,Z 

-  Lateral  model  coordinates,  in. 

xT 

-  Axial  tunnel  coordinates,  in. 

yt.zt 

-  Lateral  tunnel  coordinates,  in. 

a 

-  Angle  of  attack,  deg 

P 

-  Density,  slugs/ft3 

0c 

-  Cone  half-angle,  deg 

Introduction 

In  several  recent  tests  in  AEDC  Tunnel  B  at 
Mach  8,  there  have  been  discrepancies  observed 
between  the  measured  surface  pressure  data  and 
computed  or  theoretical  results.  For  example,  Fig. 
1 1  presents  a  comparison  of  measured  and  com¬ 
puted  pressures  on  the  Hypersonic  Lifting  Body 
(HLB)  configuration  shown  in  Fig.  2.  The  test  pro¬ 
gram  on  this  configuration  was  specifically 
designed  to  produce  an  experimental  database  for 
CFD  code  validation.  Figure  3  taken  from  Ref.  2 
represents  a  second  example  and  shows  a  com¬ 
parison  of  measured  and  computed  pressure  on 
two  axisymmetric  bodies  designed  to  produce  lin¬ 
ear  axial  pressure  distributions.  The  forebody  for 
both  of  these  axisymmetric  bodies  is  a  7-deg  cone. 


smax  ■  Circumferential  distance  around  Half  This  seemingly  consistent  bias  for  surface  pres- 
Body  sure  data  above  expected  values  was  considered 
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Surface  Static  Pressure  Ratio,  P/P< 


S/S  max  X/L 


a.  Circumferential  pressure  distribution  b.  Wind  side  pressure  distribution 

Figure  1 .  Comparison  of  computed  and  measured  pressures  on  hypersonic  lifting  body  a  =  10.1 


All  Dimensions  in  Inches 


Figure  2.  Hypersonic  Lifting  Body  configuration. 


sufficiently  troubling  to  warrant 
investigation.  This  paper 
addresses  whether,  nonuniformity 
in  the  Tunnel  B  flow  can  be 
responsible,  wholly  or  in  part,  for 
the  observed  discrepancy. 

The  general  approach  used  in 
this  study  was  to  assess  the 
effect  of  a  nonuniform  free-stream 
flow  with  computations.  This  was 
done  by  using  the  measured  tun¬ 
nel  flow  field  as  an  inflow  bound¬ 
ary  condition  for  a  CFD  solution  to 
compare  with  the  idealized  case 
where  a  uniform  free  stream  is 
assumed. 


10  15  20  25  30  35  40  M  0  5  10  15  20  25  30  35  40 

Distance  from  Nose,  X,  in.  Distance  from  Nose,  X,  in. 


a.  Surface  static  pressure  data  and  theory  for  the  b.  Surface  static  pressure  data  and  theory  for  the 
small  favorable-gradient  body  large  favorable-gradient  body 

Figure  3.  Pressure  distributions  on  axisymmetric  bodies.2 
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Two  test  configurations  were  chosen  for  study:  for  a  variety  of  configurations  and  test  conditions.  It 

a  7-deg  cone  at  zero  angle  of  attack  based  on  the  was  also  relatively  easy  to  adapt  the  code  to 

differences  shown  in  Fig.  3,  and  the  Hypersonic  accept  the  measured  free-stream  conditions. 

Lifting  Body  shown  in  Fig.  2.  Although  STUFF  has  a  nonequilibrium  capability, 

all  of  the  computations  in  this  study  assumed  per- 
Approach  fect  air. 


The  AEDC  Tunnel  B  is  a  continuous  flow 
closed-circuit,  hypersonic  wind  tunnel  with  a  50-in.- 
diam  test  section.  This  wind  tunnel,  which  is 
described  in  Ref.  3,  has  available  axisymmetric 
contoured  nozzles  to  produce  test  section  flows 
with  nominal  Mach  numbers 


In  this  paper,  it  is  necessary  to  refer  to  both  tun¬ 
nel  coordinates  and  model  coordinates.  Wind  tun¬ 
nel  coordinates  designated  by  a  subscript  T  are 
defined  in  Fig.  4a,  which  represents  a  sketch  of  the 
test  section  of  Tunnel  B.  Model  coordinate  defini- 


of  6  and  8.  A  thorough  cali¬ 
bration  of  Tunnel  B,  docu¬ 
mented  in  Ref.  4,  was  car¬ 
ried  out  in  1991.  Pitot  pres¬ 
sure  surveys,  total  tempera¬ 
ture  surveys,  and  flow  angu¬ 
larity  measurements  were 
obtained  as  a  part  of  this 
effort.  Survey  data  was 
taken  at  unit  Reynolds  num¬ 
bers  of  1.0,  2.0,  and  3.5  x 
106.  In  carrying  out  CFD 
computations,  initial  distribu¬ 
tions  of  five  flow  variables 


(Looking  Upstream) 


a.  Wind  tunnel  coordinates 


are  required:  (p,  pU,  pV,  pW, 
and  p(e+Q2/2)).  Pitot  pres¬ 
sure  and  total  temperature 
measurements  are  suffi¬ 
cient  to  determine  local  val¬ 
ues  of  p,  U,  P,  and  T.  The 
cone  (Mach-Flow  Angularity) 
probe  provides  information 
on  the  V  and  W  velocity 
components.  On  the  rake, 
there  were  42  active  pitot 
probes,  25  shielded  total 


Figure  4.  Tunnel  and  model  coordinate  systems. 


temperature  probes,  and  3  MFA  probes.  The  scar-  tions  ancj  their  relationship  to  tunnel  coordinates 
city  of  the  cone  probe  data  did  not  permit  specifica-  are  presented  in  Fig.  4b. 
tion  of  the  V  and  W  velocity  components.  These 


were  assumed  to  be  zero  at  the  inflow  computa-  The  primary  calibration  data  used  in  this  paper 


tional  plane. 


is  derived  from  the  rake.  The  42  pitot  probes  cov- 


The  CFD  solutions  for  this  study  were  gener¬ 
ated  using  the  STUFF  code  documented  in  Ref.  5. 
This  flow  solver  was  selected  because  it  has  been 
used  extensively  at  AEDC  for  supersonic  flow  com¬ 
putations  and  is  known  to  produce  reliable  results 


ered  a  range  from  ZT  =  -20  to  ZT  =  +19.8  in.  and 
the  25  temperature  probes  spanned  -16.2  <  ZT  < 
15.8.  The  rake  could  be  traversed  in  the  axial  (XT) 
direction  or  the  lateral  (YT)  direction.  Figures  5  and 
6  show  examples  of  the  pitot  pressure  and  total 
temperature  rake  data  from  a  lateral  survey  and 


3 

American  Institute  of  Aeronautics  and  Astronautics 


Yt,  in.  yt>  in- 

a.  Pitot  pressure  b.  Local  Mach  number 


Figure  5.  Rake  pitot  measurements,  Xj  =  16  in. 


sures  are  only  slightly 
affected  by  the  total  tem¬ 
perature  variations,  and 
the  computed  pressure 
results  reported  herein 
would  not  be  measurably 
changed  if  the  local  total 
temperature  values  pre¬ 
sented  in  Fig.  6  were 
replaced  by  a  constant, 
T0  =  1350”R,  the  reser¬ 
voir  temperature. 

Unfortunately,  lateral 
surveys  such  as  shown 
in  Figs  5  and  6  were  per¬ 
formed  at  only  a  few  sta¬ 
tions.  The  farthest 
upstream  station  avail¬ 
able  is  that  shown,  XT  = 
16  in.  An  axial  survey  at 
YT  =  0  extended  from  XT 


Total  Temperature  Contours,  R 


Yt,  in. 


Static  Temperature  Contours,  R 


0  10 

Yt,  in. 


=  21  to  XT  =  -14  in.  A 
difficulty  encountered 
in  carrying  out  the 
approach  was  that  for 
the  long,  slender  con¬ 
figurations  under  study, 
the  nose  of  the  model 
extended  beyond  the 
range  of  survey  data. 
Specifying  initial  condi¬ 
tions  for  the  computa¬ 
tions  was  handled  dif¬ 
ferently  for  the  two  dif¬ 
ferent  configurations. 
For  the  7-deg  cone,  the 
farthest  upstream  data 
station,  XT  =  21,  from 
the  axial  survey  was 
used  for  the  initial  con¬ 
ditions.  Moreover, 


a.  Local  total  temperature  b.  Static  temperature  since  lateral  date  was 

not  available  and  the 

Figure  6.  Rake  temperature  measurements,  XT  =  16  in.  cone  was  mounted  on 


the  Mach  number  and  static  temperature  contours 
extracted  from  that  data.  The  fact  that  the  total 
temperature  contours  are  parallel  to  the  direction  of 
rake  travel  indicates  that  these  measurement  are 
influenced  by  probe-to-probe  bias.  Numerical 
experiments  have  shown  that  model  surface  pres- 


the  tunnel  centerline, 
the  flow  field  was  assumed  to  be  axially  symmetric. 
For  the  nonsymmetric  Hypersonic  Lifting  Body,  the 
lateral  survey  data  at  XT  =  16  was  assumed  to 
apply  at  the  model  nose  at  about  XT  =  21 . 
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For  the  7-deg  cone,  computations  were  carried 
out  primarily  for  a  unit  Reynolds  number  of  1 .0  x 
106  since  much  of  the  data  of  Ref.  2  was  acquired 
at  that  value.  At  this  Reynolds  number,  viscous 
computations  can  be  performed  without  the  com¬ 
plicating  factors  of  boundary  layer  transition  and 
turbulence  modeling.  For  the  Hypersonic  Lifting 
Body  most  of  the  data  were  acquired  at  Re/L  =  3.5 
x  106/ft  with  boundary  layer  trips  installed.  The 
majority  of  the  computations  carried  out  for  this 
configuration  were  performed  assuming  an  inviscid 
fluid  but  with  initial  conditions  corresponding  to  a 
unit  Reynolds  number  of  3.5  x  106/ft. 

Figure  7  shows  measured  pitot  profiles  at  Xj  = 
21.25  for  the  lowest  Reynolds  number  and  the 
Mach  number  distributions  calculated  from  the  pitot 
pressures  and  the  measured  reservoir  pressure. 
The  drop  in  pitot  pressure  for  |  Z-j- 1  >15  is  clearly 


b.  Mach  number  distribution 
Figure  7.  Pitot  pressure  and  Mach  number  distribu¬ 
tion,  XT  =  21 .245  in.,  Re  =  1 .0  x  106  ft. 


due  to  the  tunnel  wall  boundary  layer.  The  data 
reduction  procedure,  which  assumes  inviscid  flow, 
produces  an  incorrect  increase  in  Mach  number  as 
the  wall  is  approached.  In  the  present  work,  no 
attempt  was  made  to  account  for  the  wall  boundary 
layer  in  the  computations  -  only  the  inviscid  core 
flow  was  treated.  The  test  section  area-weighted 
Mach  number  for  this  flow  condition  was  deter¬ 
mined  in  the  calibration  to  be  7.934. 

Results  and  Discussion 
Empty  Tlinnel 

Figure  8  shows  some  of  the  initial  profiles  that 
were  used  in  the  computations  for  the  7-deg  cone. 
The  flow  was  assumed  to  be  axisymmetric  and  the 
initial  values  were  interpolated  from  the  measure¬ 
ments  on  the  initial  ray  of  the  computational  grid. 
The  computational  grid  extended  from  the  tunnel 
axis  to  the  wall.  However,  as  discussed  earlier,  no 
attempt  was  made  to  model  the  tunnel  wall  bound¬ 
ary  layer.  Thus  the  data  at  what  was  judged  to  be 
the  last  point  outside  the  boundary  layer  was 
repeated  at  the  wall  for  interpolation  purposes. 

The  first  computation  was  performed  for  an 
empty  tunnel  to  compare  with  rake  calibration  data 
farther  downstream  in  the  test  section.  The  results 
of  these  computations  are  shown  in  Fig.  9.  The 
agreement  between  the  computations  and  the  data 
is  quite  good  with  the  exception  of  the  last  data 
point  shown.  The  pressure  and  density  increase  at 
this  point  is  the  result  of  a  compression  wave  that 
originates  far  upstream  in  the  nozzle  and  repre¬ 
sents  the  edge  of  the  test  rhombus.  The  good 
agreement  shown  in  Fig.  9  lends  credibility  to  the 
computational  technique  and  indicates  that  no 
unwarranted  assumptions  have  been  made. 

7-deg  Sharp  Cone 

The  next  set  of  computations  that  will  be  dis¬ 
cussed  are  for  a  7-deg  sharp  cone  positioned  on 
the  tunnel  axis  with  the  apex  located  at  XT  =  21 .25 
in.  The  computational  grid  for  this  configuration  is 
shown  in  Fig.  10  and  was  used  for  calculations 
assuming  both  an  inviscid  and  viscous  fluid.  The 
clustering  of  coordinate  lines  at  the  wall  is  neces¬ 
sary  with  a  viscous  fluid  to  resolve  the  wall  bound- 
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c.  Pressure  distribution 

Figure  8.  Initial  computational  profiles,  XT  =  21.25. 


ary  layer.  For  a  shock  capturing  flow  solver,  as 
being  used  in  the  present  study,  on  a  sharp  body, 
radial  grid  clustering  is  also  helpful  in  capturing  the 
bow  shock  near  the  apex.  To  illustrate,  Fig.  11 
shows  pressure  contours  from  a  computation  for 
an  inviscid  fluid  and  a  uniform  freestream. 


c.  Pressure  distribution 

Figure  9.  Comparison  of  computations  with  data  for 
empty  tunnel,  XT  =  -14. 

The  effect  of  the  nonuniform  free  stream  on  the 
cone  pressure  for  a  Reynolds  number  of  1 .0  x  1 06/ 
ft  is  shown  in  Fig.  12.  Figure  12a  shows  a  compari¬ 
son  of  inviscid  computations  with  the  measured  ini¬ 
tial  profiles  with  the  idealized  case  of  a  uniform  free 
stream.  The  theoretical  cone  pressure6  or  a  Mach 
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Figure  10.  Computational  grid  0C  =  7  deg. 


Figure  1 1 .  Pressure  contours  inviscid  flow,  uniform 
free  stream. 

number  of  7.934  is  also  shown  in  the  plot.  There  is 
a  definite  increase  in  cone  surface  pressure 
caused  by  the  non-uniform  free  stream.  Figure  12b 
shows  results  from  a  similar  computation  assuming 
a  laminar  boundary  layer  on  the  cone.  The  differ¬ 
ence  between  the  pressure  distributions  for  the 


a.  Inviscid  flow 

Figure  12.  Effect  of  non-uniform  free-stream  flow 
on  cone  pressure. 


b.  Laminar  flow 
Figure  12.  Concluded. 


uniform  and  nonuniform  free  stream  is  approxi¬ 
mately  the  same  for  the  viscous  and  inviscid  calcu¬ 
lations. 

Figure  13  shows  a  comparison  of  the  nonuni¬ 
form  free-stream  CFD  results  for  a  laminar  bound¬ 
ary  layer  with  a  compilation  of  data  for  a  7-deg 
cone.  The  bulk  of  the  data  shown  is  forecone  data 
previously  presented  in  Fig.  3  Some  additional 
data  at  higher  Reynolds  number  that  is  thought  to 
be  laminar  has  been  included.  The  computational 
results  represent  the  data  as  well  as  can  be 
expected. 


Figure  13.  Computational  results  with  nonuniform 
free-stream  compared  to  7-deg  cone 
data. 

The  effect  of  a  nonuniform  free-stream  on  the 
inviscid  pressure  distribution  at  a  unit  Reynolds 
number  of  3.5  x  1 06/ft  is  shown  in  Fig.  14.  Although 
the  computations  assume  an  inviscid  fluid,  the  ini¬ 
tial  flow  distribution  depends  on  tunnel  Reynolds 
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0.05 


5  10  15  20  25  30 

X,  in. 

a.  Uniform  initial  profiles 


o. 

o. 


0. 

Q. 

o 

0. 

0. 


X,  in. 

b.  Measured  initial  profiles 
Figure  14.  Effect  of  initial  profiles  on  cone  pressure 
at  Re  =  3.5  x  106ft. 

number.  This  computation  uses  initial  data  derived 
from  Figs.  5  and  6  and  takes  into  account  the 
three-dimensional  character  of  the  flow.  The  pres¬ 
sure  distribution  for  the  measured  initial  profiles  is 
virtually  identical  to  that  for  a  uniform  free  stream. 
Thus  the  effect  of  a  nonuniform  free  stream 
appears  to  be  much  less  at  the  high  Reynolds 
number  than  at  low  Reynolds  number.  Presum¬ 
ably,  the  smaller  displacement  thickness  on  the 
nozzle  wall  at  high  Reynolds  number  results  in 
somewhat  better  flow  quality. 

Hypersonic  Lifting  Body 

Results  were  computed  for  the  HLB  configura¬ 
tion  for  the  range  of  positive  angles  of  attack  cov¬ 
ered  by  the  experimental  program.  Three  cross- 
sectional  planes  from  an  example  grid  used  for  the 
computations  at  an  angle  of  attack  of  10  degs  are 
shown  in  Fig.  15.  The  outer  boundary  of  the  grid  is 


a  cylinder  of  radius  15  in.  that  is  parallel  to  the  tun¬ 
nel  axis.  A  30-in.-diam  circle  includes  the  region 
covered  by  the  rake  survey,  except  for  the  areas 
where  the  tunnel  wall  boundary  layer  intrudes. 


a.  Cross-sectional  grid  -  X  =  0.6893  in. 


Y,  in. 

b.  Cross-sectional  grid  -  X  =  18.09  in. 


c.  Cross-sectional  grid  -  X  =  40.0  in. 
Figure  15.  Example  cross-sectional  grids  for  HLB  at 
a  =  10  deg. 
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Constructing  a  grid  for  the  model  in  the  tunnel 
required  knowledge  of  the  location  of  the  center  of 
rotation  (see  Fig.  4)  in  model  coordinates  and  tun¬ 
nel  coordinates.  The  test  project  engineer  provided 
these  values  (Xc  =  34,  Yc  =  0,  Zc  =  0,  XTC  =  -18, 
Ytc  =  0,  ZTC  =  0). 

Figure  16  shows  inviscid  computational  results 
for  both  a  uniform  and  nonuniform  free  stream 


c.  X/L  =  0.80 

Figure  16.  Circumferential  pressure  distributions, 
inviscid  computations,  a  =  10  deg. 


compared  with  measured  pressure  data  for  10-deg 
angle  of  attack.  There  is  very  little  difference 
between  the  two  computations,  both  being  about 
10  percent  below  the  data  on  the  wind  side.  Similar 
computational  results  for  a  1 5-deg  angle  of  attack 
are  presented  in  Fig.  17.  For  this  case  there  is  a 
noticeable  increase  in  wind-side  pressure  due  to 
the  nonuniform  free  stream.  This  increase  is 


b.  X/L  =  0.70 


c.  X/L  =  0.80 

Figure  17.  Circumferential  pressure  distributions, 
inviscid  computations,  a  =  15  deg. 
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thought  to  be  due  to  the  fact  that  the  nose  and  fore¬ 
body  of  the  vehicle  at  this  angle  of  attack  were 
positioned  far  enough  from  the  tunnel  centerline  to 
encounter  a  slightly  higher  Mach  number  flow  (cf. 
Fig.  5).  The  computed  values  on  the  wind  side  are 
still  significantly  below  the  data. 

Unfortunately,  the  computations  could  not  be 
completed  over  the  entire  vehicle  for  20-deg  angle 
of  attack  due  to  negative  densities  encountered  on 
the  lee  surface.  Comparisons  of  inviscid  computa¬ 
tions  with  pressure  data  at  X/L  =  0.25  and  X/L  = 
0.55  shown  in  Fig.  18  reveal  a  noticeable  effect  of 
the  wind  side  due  to  the  nonuniform  free  stream. 
The  computational  results  for  a  nonuniform  free 
stream  are  definitely  approaching  the  data  as 
angle  of  attack  increases.  This  is  almost  certainly 
due  to  the  radial  distributions  of  Mach  number 
shown  in  Fig.  5. 


b.  X/L  =  0.55 

Figure  18.  Circumferential  pressure  distributions, 
inviscid  computations,  a  =  20  deg. 


Summary  and  Conculsion 

The  results  of  this  study  indicate  that  flow  non¬ 
uniformities  are  at  least  partially  responsible  for 
discrepancies  observed  between  measured  and 
computed  surface  pressures  at  low  Reynolds  num¬ 
ber  (Re/L  =  1.0  x  106/ft)  on  slender  axisymmetric 
bodies  at  zero  angle  of  attack.  There  is  no  corre¬ 
sponding  effect  at  high  tunnel  Reynolds  number 
(Re/L  =  3.5  x  106/ft).  Presumably,  the  flow  quality 
is  better  at  high  Reynolds  number  due  to  a  thinner 
tunnel  wall  boundary  layer.  Even  a  high  Reynolds 
number,  however,  for  angles  of  attack  that  displace 
the  nose  of  the  vehicle  more  that  6-7  in.  from  the 
tunnel  centerline,  there  is  a  small  increase  in  sur¬ 
face  pressure  caused  by  a  radial  variation  in  Mach 
number.  Presumably,  yaw  or  lateral  displacements 
would  cause  a  similar  effect. 

The  effect  due  to  radial  Mach  number  variations 
is  not  sufficient  to  explain  fully  the  pressure  dis¬ 
crepancies  HLB  observed,  although  the  trend  with 
increasing  angle  of  attack  is  in  the  right  direction. 
Other  possible  sources  of  the  discrepancy, 
although  unlikely,  are  (1)  model/sting  deflections, 
(2)  flow  angularity  (neglected  in  this  study),  (3)  dif¬ 
ferences  between  the  computational  model  and 
the  tunnel  model  geometry,  (4)  a  consistent  bias  in 
the  pressure  measurements,  and  (5)  errors  in  the 
computations.  One  is  reluctant  to  focus  on  compu¬ 
tational  errors  because  other  investigators  using 
different  codes  have  found  similar  results. 
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